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PRINCIPLES OF MOMENT DISTRIBUTION APPLIED TO STABILITY OF STRUCTURES 

COMPOSED OF BARS OR PLATES 

By Eugene E. Lundqtjist, Elbkedge Z. Stowell, and Eyan H. Schefettb 


SUMMARY 

The principles of the Gross method of moment distribution , 
which hare previously been applied to the stability of structures 
composed of bars under axial load, are applied to the stability 
of structures composed of long plates under longitudinal load. 
A brief theoretical treatment of the subject, as applied to struc- 
tures composed of either bars or plates, is included, together 
with an illustrative example for each of these two types of 
structure. An appendix presents the derivation of the formulas 
for the various stiffnesses and carry-over factors used in solving 
problems in the stability of structures composed of long plates. 

INTRODUCTION 

The usual procedures for calculating critical buckling 
loads for the members of complex structures are often some- 
what involved and are not easily reduced to a set of routine 
calculations. Many practical engineers, as a consequence, 
do not attempt to calculate critical buckling loads. 

One approach to the solution of problems in the stability 
of structural members that is purely engineering in character 
and that lends itself to simplified calculations is provided by 
use of the principles of the Cross method of moment dis- 
tribution (reference 1). The theory of moment distribution, 
originally devised as a rapid method of stress analysis, 
describes how the resistance to an external moment, applied 
at any joint in a structure composed of bars, is distributed 
throughout the structure in accordance with the resistance 
of the various joints to rotation. The original theory of 
Cross was modified by James (reference 2) to take into 
account the possibility of axial load in the members. 

The modified theory of James has already been applied 
in reference 3 to the study of the stability of structures 
composed of bars under axial load. Because of the funda- 
mental character of the quantities used in the method of 
moment distribution and of the formulas associated with 
them, it is possible by suitable definition of the quantities to 
apply an analysis exactly like that of reference 3 to the study 
of the stability of structures composed of plates under 
longitudinal load. 

The present report gives a generalized derivation of the 
formulas, applicable to both bar and plate structures. The 
evaluation of various quantities for structures composed of 
bars was given in reference £. The corresponding evaluation 


of the quantities for structures composed of plates is given 
in an appendix to this report. 
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modulus of elasticity 
load on structure 
rotation of joint 
deflection 

series stability factor 
modified stiffness stability factor 


BARS 


effective modulus of elasticity for stresses beyond the 
elastic range 

moment of inertia of cross section about an axis per- 
pendicular to plane of bending 
area of cross section 


(VI) 


radius of gyration 

length of bar 
axial load in bar (absolute value) 


fixity coefficient in column formula 


stiffness factor 





PLATES 


effective plate modulus for stresses beyond the elastic 
range 

Poisson’s ratio 

half wave length of buckles in longitudinal direction 
width of plate 
thickness of plate 

(Py \ 

1 2 ( 1 — u 1 ) / 

effective flexural stiffness of plate for stresses beyond 
the elastic range 
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a longitudinal compressive stress in plate 

k=-^-L<r (always positive) . 

irD 

M bending moment 

M 0 amplitude of sinusoidally distributed moment 
e restraint coefficient 
w deflection normal to plane of plate 

SUBSCRIPTS 

i initial value 
cr critical 
eff effective 
F flange 
\V web 

DEFINITIONS 

Member. — The word “member” is used in this report to 
indicate either a bar or an infinitely long, flat, rectangular 
plate. 

Joint. — A joint in a structure composed of plates, by 
analogy to a joint in a structure of bars, is defined as the 
entire length of the intersection line between two or more 
joined plates. 

Stiffness and oarry-over factor. — If a bar is on unyielding 
supports at each end, the moment at one end necessary to 
produce a rotation of one-fourth radian at that end is called 
the stiffness of the bar and the ratio of the moment developed 
at the far end to the moment applied at the near end is 
called the carry-over factor of the bar. 

In order to write similar definitions of stiffness and carry- 
over factor for plates, it is necessary to include a statement 
showing how the moment is distributed along tho edges of 
the plate. The solution of the differential equation for the 
critical compressive stress of an infinitely long plate with 
given edge restraints reveals that, when the plate buckles, 
the moments and the rotations at both edges of the plate 
vary sinusoidally along the edges and are in phase with each 
other. The ratio of moment per unit length at any point 
along the edge to the rotation at that point is therefore 
constant along the edge for a given wave length. The 
following definitions of stiffness and carry-over factor for 
plates may therefore be written: 

Stiffness — If an infinitely long flat plate, is under 
longitudinal compression with one unloaded edge on 
an unyielding support, the ratio of moment per unit 
length at any point along this unloaded edge to the 
rotation in quarter-radians at that point when the 
moment is distributed sinusoidally is called the stiffness 
of the plate. 

Carry-over factor— The ratio of the moment per 
unit length developed at any point along the far un- 
loaded edge to the applied moment per unit length at 
the corresponding position along the near unloaded 
edge is called the carry-over factor of the plate. 

The foregoing definitions make it possible to use various 
stiffnesses and carry-over factors in a similar manner for 
both bars and plates. 

The symbols used to designate the stiffness and carry-over 
factor for the different types of support and restraint at the 
far end or edge are given in the following tabl,e: 


Stiffness 

Carryover 

factor 

Conditions at far and or edge 

8 

S 1 s 
.£11 
£311 

£3V 

C 

Ci 

Cii-0 

Cui-o 

C"— 1 

Far end or edge supported and filed against rotation. 

Far end or edge supjwrtod and elastically restrained 
against rotation. 

Far end or edge supported with no restraint against 
rotation. 

Far edge free (no support and no restraint against rota, 
tion). This condition is not used in connection with 
bars. 

Far end or edge supported and subjected to moment 
equal and opposite to that applied at near end or edge. 


The quantities S l , C 1 , S n , C’ 11 of this paper correspond to 
S', C , S" , C", respectively, of reference 3. 

The stiffness of a bar computed according to the definition 
used herein is one-fourth that computed according to the 
definition used by Cross (reference 1). In moment distribu- 
tion the relative, not the absolute, values of stiffnesses of 
the members are of importance. The foregoing definition 
was selected so that the stiffness of a bar of constant cross 
section with no axial load and fixed at tho far end would be 
ElfL instead of 4 EIIL. 

Sign convention. — A clockwise moment acting on the end 
of a bar or at any station along the side edge of a plate is 
positive and causes positive rotation at that end or station. 
An external moment applied at a joint is considered to act 
on the joint; a counterclockwise moment acting on a joint 
is positive. 

CRITERION FOR STABILITY 


It is assumed that all members in a structure composed of 
bars lie in the plane in which buckling occurs and that the 
joints of the structure are held rigidly in space but are free 
to rotate subject to the elastic restraint of tho connecting 
members. Similarly, in a structure composed of plates, it is 
assumed that the joints between plates, or between plates 
and longitudinal restraining members, remain in their 
original straight lines but are free to rotate subject to the 
elastic restraint of the connecting members. 

In the discussion that follows, either of two criterions for 
stability may be used. For each criterion, the stiffness and 
carry-over factor arfe functions of the axial load in the bar 
or the longitudinal load in the plate. (Sec references 2, 3, 
4, and 5.) 

Stiffness criterion for stability. — From a structure of many 
members the section comprising one joint shown in figure 1 
is considered. Figure 1 may be interpreted as being either 
a plan view of a structure composed of bam or an end view 
of a structure composed of long plates. An external moment 
of —1 is assumed to be applied at the joint t. If the struc- 
ture is composed of plates, this moment is the external 
moment per unit length at the station under consideration. 
Because_the angles between members at the joint aro pre- 
served and the rotations of all members at the joint must 
therefore be equal, the moment of 1 added to balance this 
joint is distributed among the members in proportion to their 
stiffnesses, as follows: 


5 % 

2SV 


to member iji 


SS 1 ,, 


to member ij t 
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and so forth. The moment-distribution analysis is now 
complete as far as moments at joint rare concerned. 

For stability, the moment in the members must be finite. 
The stiffness criterion for stability is therefore 


( 1 ) 

The condition of neutral stability gives the critical buckling 
load for the structure and is obtained by setting. the stiffness 
stability factor equal to zero, or 

2S I y =0 (2) 

In the general case there is more than one critical buckling 
load; thus, satisfaction of equation (2) is insufficient for the 
solution of a given stability problem. Instead, the lowest 
load that satisfies equation (2) must be calculated and 
compared with the load for which the structure is designed. 
Only if this low'est critical load is greater than the design 
load is the structure stable. 


h 



According to the definition of stiffness, the moment dis- 
tributed to any member must be the rotation of the joint 
multiplied by the stiffness of the member. Hence 9, the 
rotation expressed in quarter-radians of joint i caused by 
the external moment — 1 , is 

0=v*i (3) 

ij 

Equation (3) will be used under the section “Method of 
Making Preliminary Estimate of the Critical Load.” 

Series criterion for stability. — In a structure of many 
members, the section comprising two joints shown in figure 2 
is considered. An external moment of —1 is assumed to 
be applied at joint ?. If the structure is composed of plates, 
this moment is the external moment per unit length at the 
station under consideration. By a moment-distribution 
analysis of reference 3, the total moment in members ih at 
joint i is 

y.Q i 

S«+2S I u (1+r+rl+r * + ' * 0 


or 

2 S^jjt 1 

Sij~r 2S 1 « 1 — r 

where 

SjjCjj SjjCji 

Su+ZS'u Sv+ZlS 1 jt 


(4) 

(5) 


For stability, the total moment in members ih must be 
finite. The series criterion for stability is therefore 


r<l (6) 

The condition of neutral stability gives the critical buckling 
load for the structure and is obtained b;y setting 

r=l (7) 

The same considerations that apply to the stiffness 
criterion for stability also apply to the series criterion for 
stability. The lowest load that satisfies the equation for 
neutral stability (in this case, equation (7)) must be cal- 
culated and compared with the load for which the structure 
is designed. If this lowest critical load is greater than the 
design load, the struct ureis stable. 

According to the definition of stiffness, the total moment 
in members ih at joint i must be the rotation of joint i 
multiplied by the total stiffness of members ih. Hence 6, the 
rotation in quarter-radians of joint i caused by the external 
moment —1, is: 

1 — r (8) 


Formulas (2) and (7) are both derived in reference 3. 
Whether formula (2) or formula (7) is to be used will depend 
upon the particular problem. In cases in which the structure 
is symmetrical about a joint, the expressions concerned with 
the stiffness criterion usually involve fewer calculations; 
when the structure is symmetrical about a member, the 
formulas concerned with the series criterion offer certain 
advantages. 

Stiffness criterion, for stability when structure is sym- 
metrical about a member. — A modification of the stiffness 
criterion in which the values of S lv are used is sometimes 
convenient when the structure is symmetrical about a 
member, as shown in figure 3. When this criterion is used, 
opposing unit moments are applied at the two ends or edges 
of the member about which the structure is symmetrical. 
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The stiffness stability factor of equation (2) for the joint ! 
in figure 3 is then written: 

+ 25^=0 (9) 

An illustration of the use of this special application of the 
stiffness criterion in a plate problem is included in the 
section “Examples.” 


CARRY-OVER FACTOR AND STIFFNESS 


In order to check the stability of a group of structural 
members by use of the equations previously given, additional 
equations for the carry-over factor and stiffness are required. 

The member ij shown in figure 4, on an unyielding sup- 
port at i and elastically restrained at j by members jk is 
considered. The members jk are also elastically restrained 
at their far ends k. By a momont-distribution analysis 
(reference 3) it follows that the carry-over factor C l tj is 


rt i n 

and the stiffness S l (} is 

qi _ 


( 10 ) 


(ii) 


Substitution of equation (10) in equation (ll)jgives 


S\ } => 


S »„ 

i ‘ n n __ 

CliCi >S n ]t +xS I Jk 


( 12 ) 


For member ij, the limiting values of the carry-over 
factor and of stiffness given by equations (10) and (12), 
respectively, are obtained as follows: When the far end j is 
pinned, there is no elastic restraint at j and 25' 1 Jt =0. For 
this limiting condition, C 1 tJ =C ll tj =0, and S 1 il =S 1I ii . When 
the far end j is fixed, there is complete restraint at j and 
2(S I ^=co . For this limiting condition, CP if =C (1 and S 1 ij =S t} 
where 


Stl ~i 


(13) 


A similar equation, which expresses S lv tJ in terms of S 11 ^ 
and C }t , can be obtained from equation (11) as follows: If the 
restraint at the far end is such that C l i} =— 1, there must be, 
at the far end, a moment of the same magnitude but opposite 


in direction to that applied at the near end. If, therefore, 
(Pt) in equation (11) equals —1, S 1 ^ becomes S xy tj, where 

<«> 

The expressions used for the computation of numerical 
values of S, C, S 1I , <S IH , and S lv for plates are given in the 
appendix. 

Up to this point, all the equations in this report are general. 
In nearly all cases encountered in practice, however, the 
cross section and axial load do not vary along the length of 
each member. For this special case, C ti =Cjt, S t j—Sjt, and 
so forth. In practical problems the numerical values for 
these quantities are obtained by use of tables. Such tables 
are given for bars in reference 4, where the argument is 
( Ljj ) erf , and for plates in reference 5, where the arguments 
are k and X/6. 





METHOD OF MAKING PRELIMINARY ESTIMATE OF THE 
CRITICAL LOAD 

In order to determine the lowest critical load for tho 
structure, it is necessary to test either equation (2) or equa- 
tion (7) for neutral stability for different assumed loads. 
The lowest load that satisfies either equation is the critical 
load for the structure. If evaluation of the stiffness or the 
series stability factors has required lengthy compulations and 
if all the assumed loads for which these factors have been 
evaluated arc less than the lowest critical load, as evidenced 
by the fact that 2S 1 ,, remains positive or that r remains less 
than unity, a method utilizing tho work already clone may 
be used to estimate the critical load. This estimated load may 
then be used as a trial load in equation (2) or equation (7). 

The method of estimating the lowest critical load is based 
upon principles used in the analysis of experimental obaer- 
vationsJn problems of elastic stability (references 6 and 7). 
Southwell (reference 6) mentioned that the unavoidable im- 
perfections in practical structures prevent the realization of 
the concept of a critical load at which deflections begin. 
Instead, the initial deflections present in practical structures 
steadily grow with increase in load and, according to the 
usual theory, the deflections become infinite as tlu* critical 
load is approached. 
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The general relation between load and deflection for 
problems of elastic stability (reference 7) shows that if 
(y—yi)l(P—Pi) is plotted as ordinate against y—y L as 
abscissa, the curve obtained when P approaches P„ is essen- 
tially a straight line of which the inverse slope is P„—P x . 
Here y is the deflection at load P in a member, yi and Pi 
are initial values of y and P, respectively, P„ is the lowest 
critical load, and 

Pi<P<P er 

If simultaneous readings of load and deflection recorded in 
a test are plotted as described with any load P as the initial 
reading, the value of Pa-— Pi is readily obtained. The value 
of P„ is then given by the relation 


Pcr=(P c -Pi)+Pi (15) 

The relation between load and deflection can also be 
applied to load and rotation of a joint provided that there is 
an initial rotation of the joints. The initial rotation is 
obtained by the application of the external moment —1 at 
some joint, after which the load on the structure is applied. 
As the lowest critical load is approached, the rotations 
become infinite. 

If the distribution of the loads throughout the structure 
does not change as the total load IT increases, the axial or- 
longitudinal load in each member is proportional to IT. 
'If (0— IVi) is plotted as ordinate against 6—di as 
abscissa, the curve obtained when IT approaches W a is 
essentially a straight line with inverse slope TT„— TTi, 
where 0 is the rotation of a joint under the external moment 
— 1 when load IT is on the structure, 0i and 11^ are initial 
values of 0 and 11', respectively, !T cr is the lowest critical 
load, and 


H I<1H < Cit'w 

When simultaneous values of load and rotation are plotted 
as described with IFi as the initial load, the value of IFi 
is easily obtained. The value of W„ is then given by the 
equation 

ll' cr = (!!'„- H'O + II'i (16) 


The procedure to be used in estimating the critical load 
for a group of structural members is as follows: 

1. For each of the loads IT assumed in the application of 
one of the stability criterions (equation (2) or equation (7)) 
to a joint, calculate the rotation 0 of this joint by means of 
equation (3) or equation (8). 

2. Designate the lowest assumed value of IT and the 
corresponding value of 0 as IFi and 0 If respectively. 

3. Plot the curve of (0— 0i)/(TT— H'O as ordinate against 
0—0i as abscissa and estimate W er from the slope of the 
resulting line. If the curve obtained is not essentially a 
straight line, successively higher values of the assumed loads 


W should be designated TTi and the value of re-estimated. 
The accuracy of the estimated value of TTo- is improved as 
both TT and Tt'i approach TT„. 

An example of the application of this method for predict- 
ing the lowest critical load is given in reference 8. 

As applied to a structure of plates, this method gives a 
critical load for some particular value of the half wave 
length X. The value of TT CT that satisfies equation (2) or 
equation (7) and is a minimum with respect to X must 
finally be found as in the example, given subsequently 
herein, in which the use of this method of estimating the 
critical load for a given wave length was not required. 

DISCUSSION OF METHODS 

Each of the two equations for neutral stability contains 
the stiffness of certain members elastically restrained at their 
far ends or edges by other members. These other members 
may also be elastically restrained at their far ends or edges 
by still other members, and so on. By successive applica- 
tion of equation (12) the restraining effect of all the members 
in the structure can be considered. 

In practical calculations for structures composed of bars, 
modification of the actual structure by the introduction of 
pins at certain joints is usually necessary. It has sometimes 
been the custom to consider only one member elastically 
restrained at the ends by the adjacent members, which are 
assumed to be pinned at the far ends. The calculation of 
W„ by use of small groups of members in this manner is 
quite inadequate. Treatment of much larger groups of 
members in one calculation is necessary if a reasonably 
accurate value of TT^ is to be obtained. 

If the stresses in any of the members of a structure are 
beyond the elastic range, the reduction of the modulus of 
elasticity at these stresses must also be considered. Dis- 
cussions of this reduced modulus for structures composed of 
bars are given in references 3 and 8. References 9 and 10 
discuss the reduced modulus of elasticity for plates at high 
stresses. 

EXAMPLES 

Structure composed of bars. — The example of a structure 
composed of bars presented herein is identical with that 
given in reference 3; for the solution of the problem, the 
tables of reference 3, rather than the more extensive tables 
of reference 4, were used. 

A continuous member of 1026 steel is to be designed to 
carry the loads shown in figure 5. For simplicity, the same 
cross section will be used in all spans. 


£ 

V 


AxicH load in pounds: T, tension ; C, cornpressian 
O 9940 C 86(0 T 9940 C 86/0 T 9940 C 0 

A A A A A A A 

z a b c d e f 


f — GO 


5 at 50*250 


60 — > 


Figcei e.— Illustrative bar problem. 
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The usual column formulas for 1025-steel tubes are: 


For ^<124, 

•^- 36 , 000 - 1.172 1 (jj 

For ->124, 

P 

P 276X10* 



(17) 


( 18 ) 


2. For each assumed axial load in the selected member, 
calculate the corresponding axial load in every other member. 
In this problem the axial load in all compression members 
is the same and the axial load in the tension members is 
0.866 times the axial load in the compression members. 

_3. For each load in each of the members, calculate P/A, 
E, and ( L/j) tff . In this problem, E is obtained from equation 
(17) by methods outlined in reference 3, or 



It is desired that Ljp bo less than 124. Equation (17) there- 
fore is used and, on the assumption that c=2, a tube of the 
following dimensions is selected as a trial design for com- 
pression members za, be, and de. 


Diameter, d in... 1. 625 

Wall thickness, t in. 0. 065 

Area, A _sq in... 0. 3186 

Moment of inertia, I in. 4 .. 0. 09707 

According to the problem, this tube is used as a continuous 
member from y to / (fig. 5). 

In order to check the stability of the tube selected in the 
trial design, the critical buckling load will be calculated and 
compared with the loads given in figure 5. The axial load 
in the tension spans is assumed to be always 8610/9940 or 
0.866 times the axial load in the compression spans. This 
assumption conforms to the condition that the forces in all 
members increase in the same ratio as the load on the 
structure increases. 

Both the dimensions and the loading of the member shown 
in figure 5 are symmetrical about span be. It is therefore 
convenient to determine the critical buckling load by use of 
the series criterion for stability. If the unit external moment 
to be applied is at joint b, the series stability factor is given 
by equation (5) .with the summation signs omitted. If the 
symmetry about span be is considered, the series stability 
factor becomes 


where 


j 

r 7(’S»c+s i «,) 2 




s 11 , 


at 


s 1 * 


1 pi & 

1 C ** S"<~+S* cf 


( 19 ) 


In the equation for S 1 <u it is assumed that the ends at y and/ 
are pinned. 

The detailed procedure of calculating the critical buckling 
load is as follows: 

1. Assume a series of values for the axial load in one of 
the members. In order that the values of load be reasonable, 
a compression member should always be selected and the 
values of the axial load for this member computed from the 
column formula by use of a series of values of c. In this 
problem, compression member be is selected and the column 
formula is equation (17). 


4. For each load in each of the members, determine the 
value of the terms required to evaluate equation (19), by use 
of the tables of reference 3 or 4. 

5. The assumed load that gives r= 1 is the critical buckling 
load. 

The results of this procedure as applied to the problem of 
figure 5 are given in table I ; the values of c in the fust column 
are given for reference only and, as stated in paragraph 1 of 
the foregoing procedure, were so assumed that a scries of 
reasonable values for the axial load P in the compression 
member be could be obtained. In the last column of table I 
are given the values of r corresponding to the assumed 
values of c. As the value of c increases from 1.4 to 2.6, (he . 
value of r increases from 0.133 to 1.63. If the data of table I 
are plotted, it is found that when r— 1 the lowest critical 
buckling loads for the trial design are 

za, be, and de... 10,260 compression 

ab and cd 8,890 tension 

These critical loads are greater than the loads to which the 
respective members are subjected (see fig. 5). The tube 
selected for the trial design is therefore stable and the margin 
of safety for the system is 


10260 8890 

9940 8610 


1 = 0.03 


A single margin of safety is obtained for the whole system 
regardless of which member is used for its calculation because, 
when the critical load is reached, all members deflect. 

More than one type of instability is possible, theoretically; 
therefore, as the loads P increase, there is more than one 
value of P for which r— 1. ‘ (See table I.) For each type of 
instability there is a corresponding critical load. In design, 
however, the lowest critical load should be calculated and 
compared with the loads given in the problem. 

Table I shows further that, for values of c between 1 .4 and 
1.5, the value of S changes from positive to negative. 
According to the stiffness criterion for stability, this change 
of sign means that members de and ef, considered alone, have 
changed from stable to unstable. It is also noted that S 1 # 
changes from positive to negative for values of c between 2.0 
and 2.7 ; members cd, de, and ef, considered alone, have there- 
fore changed from stable to unstable, but at a much higher 
load. The change from stable to unstable for all members 
occurs for values of c between 2.5 and 2.6 when r= 1. 
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Structure composed of plates. — The critical compressive 
stress for local instability of a 24S-T aluminum-alloy 
Z-section column with the cross-sectional dimensions shown 
in figure 6 is to be determined. 

It is convenient in symmetrical plate problems of this 
type to use the modification of the stiffness criterion for 
stability previously discussed. If opposing unit external 
moments are applied at the joints between the web and 
the flanges, the stiffness stability criterion, as given by 
equation (9), is 

ZS l v =S m r +EF v w =0 (20) 

where the subscripts F and IT refer to the flange and the 
web, respectively. 



The tables of reference 5 give the values of S in and S iy 
in the dimensionless form S lu [(Dlb) and SP^KD/b) rather 
than directly. It is therefore desirable to write equation 
(20) in the form 



If this equation is divided by 



it becomes 



Because T) w lDr=(tw/tr) a , the stability criterion may be 
written in terms of the modified stiffness stability factor U, 
as 



The detailed procedure of calculating the critical com- 
pressive stress is as follows: 

1. Compute the ratios t w jt F , b F jb Wl and bw/t w . 

2. Assume a value of \/b F . 

3. Compute X/S^ from the equation 

bw bp bw 

4. Assume a series of values of k F and, for each value of 
k F , compute k w from the equation 



The indicated procedure is adopted as being somewhat more 
convenient than assumption of the stress and computation 
of the corresponding values of k F and k w . It is permissible 
to compute from the given equation even though the 
stress is beyond the elastic range, because the stress and^ 
thus, by assumption, the effective plate modulus are the 
same in the web and the flange. 

5. Evaluate the modified stiffness stability factor U 
from equation (21) and the tables of reference 5. 

6. Plot U against k p or k w and note the value of k for 
which U is equal to zero. 

7. Repeat steps 2 to 6, assuming different values of \/b F . 

8. Plot values of k F for (7=0 against \/b F (or k w for 
(7=0 against \/b w ) to determine the minimum value of 
k F (or k w ). 

9. With this minimum value of k, evaluate the critical 
stress from the formula (see definition of k), 


" ¥T 


which may be written, for the web, 


„ k Et it 1 

" 12 (1-/)^ 

or, for the flange, _ 

kpTpEtp 2 

* 12 ( 1 — fP)b F * 


( 22 ) 


(23) 


The value of <t„ will be the same regardless of whether 
equation (22) or (23) is used. 

The results of this procedure as applied to the problem 
of figure 6 are given in table II. The values of k w for (7=0 
in the last column of table II were determined according 
to step 6. If these values of k w are plotted against \/b w 
(step 8), the minimum value of k w is found to be about 2.9. 
(See fig. 7.) The critical compressive stress for local buckling 
of the section shown in figure 6 is then, from equation (22) , 

2.9 X 9.87 X 10.6 X10 7 , . , 

Vcr= — 12X0 91 X (40) 2 — =17 > 400 pounds per square inch 


This method provides a relatively simple means of predict- 
ing the critical-stress values for columns of Z-section and 
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Figure 7.— Plot of kw against \Jbr for plate problem. 

other simple cross sections, such as I-, channel, and rectan- 
gular-tube sections. Charts giving the values of k deter- 
mined by this method which were prepared for wide ranges 
of the dimension ratios are presented in reference 11 for 
columns of I-, Z-, channel, and rectangular-tube section. 


An alternate method of solution for problems of this type 
makes use of the charts of references 9 and 10 and the tables 
of reference 5. An assumption is made as to whether the 
flange or the web will be primarily responsible for instability. 
If the flange is expected to be primarily responsible, the value 
of S lv w for the web is determined from the tables of reference 5. 
This value is then used in computing the restraint coeffi- 
cient e (reference 9 or 10), and the value of k is found from 
figure 3 of reference 9. Because it is necessary to assume a 
value of k and A/6 in order to determine S IV W , the method 
will obviously involve a trial-and-error procedure. Further- 
more, if repeated calculations show that S 1V is negative, the 
assumption that the flange would bo primarily responsible 
for instability is incorrect. In this case, it will bo necessary to 
evaluate S ni F and to determine it from figure 3 of reference 10. 
A detailed example of the application of this method is given 
in reference 11. 

Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronauticb, 
Langley Field, Va., July 15, 1948. 




APPENDIX 

DERIVATION OP STIFFNESSES AND CARRY-OVER FACTORS 


PLATE UNDER COMPRESSION 


In order to apply the method of moment distribution in 
any form, the values of stiffnesses and carry-over factors are 
required for the members in question. Formulas for the 
evaluation of these quantities for bars were developed in 
reference 3. This appendix gives the corresponding deriva- 
tion of the formulas for plates; the sign convention used, as 
distinguished from that given in the section on “Definitions,” 
corresponds to that of reference 12, in which deflections w are 
positive downward and a moment is positive if it produces 
compression in the upper fibers. 

General deflection surface of a plate buckled under com- 
pression. — Before the values of stiffness and carry-over factor 
for flat plates under various conditions of edge restraint may 
be computed, the deflection surface of a flat plate buckled 
under a compressive load with a moment applied along one 
unloaded edge must be described. 

An infinitely long flat plate under longitudinal compression 
is shown in figure 8 with coordinate axes. For equilibrium of 
an infinitesimal element of the plate, the following equation 
must be satisfied (reference 12, p. 324): 


dhp , 0 dho r dhu , bhv 
dr 4 6 2 dx*~ 


(Al) 


On the assumption that the plate is infinitely long in the 
direction of x, the conditions at the ends do not matter; the 
solution of equation (Al) is therefore taken in the form 


W=f(y) cos y 


(A2) 


The unknown function f(y) may be determined by sub- 
stituting the expression for w into equation (Al). It is 
found that the function/ must satisfy the equation 

f? + (£~w) /=0 ^ 

Equation (A3) is an ordinary differential equation of the 
fourth order, the solution of which is 

f—Ci cosh °j+(h sinh y+c 3 cos y-f-c 4 sin (A4) 

where Ci, Cj, c s , and c* are arbitrary constants and 



+x 



WW-rt-V* 

The deflection surface of the plate is now found by ^sub- 
stituting this result for / in equation (A2): 

w=(ci cosh ^y-f-e 2 sinh ^+e 3 cos ^+c 4 sin cos ^ (A5) 

In this solution, four conditions may be imposed along the 
unloaded edges to fix the four constants. One of the four 
conditions will always specify the presence of a moment 

M 0 cos along the near edge, and another will specify that 

the deflection w along this edge is zero. The remaining two 
conditions will be varied to suit the conditions at the far 
edge of the plate of which the stiffness is being computed. 
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Stiffness of a plate with far edge fixed. — Figure 9 shows 
a flat rectangular plate under compression with a moment M 

applied along one edge at y= — ^ and with complete restraint 

against rotation along the parallel edge at y =^- The stiffness S 
of the plate is defined as 




CA6) 


where (0) 4 is the rotation of the edge at y— — w expressed 

f—T 

inquarter-radians. 

The general expression (A5) for the deflection of the plate 
must be specialized to the case of figure 9 in which the 
boundary conditions are: 


w ,-4 =0 - 


-Kv+'I?) ,- M=M ' c ° s ” 


r— 


w= 


M 0 b 2 


(ttr 0 

(a tanh ^+0 tan 


(A7) 

(AS) 

(A9) 
sinh ^ sin ^ 



After determination of the arbitrary constants in equa- 
tion (A5) by use of these boundary conditions, the deflection 
surface for the case of figure 9 is found to bo 


. , a . 0 
sinh ^ sm ^ 


D(c?+1 3 2 ) 




(a cotli ^ — 0 cot ^ 


cosh 


ocy 


cos - 1 - 


§y'S 


cosh ^ cos ~ 


oc tanh ^+0 tan |+« coth 0 cot | a tanh 0 tan |+a coth 0 cot ^ 


cos 


TTX 


(A10) 


From this deflection surface there is obtained 


( 0 ) 


=4<^ 


Mb 


r— -I | Z>(a*+i 3 s ) 


(All) 


a'tanh ^ + 0 tan § a coth ^ — 0 cot ^ 


where 0 is'expressed in quarter-radians. Substitution in equation (A6) gives 

r 


-4? 


!/+(!)’: 


i 


a , a . 0 0 a . a 0 0 

^2 tanh 2+2 tan | ^coth^-^cot^ 


S 11 


(A12) 


Carry-over faotor of a plate with, far edge fixed. — The 
carry-over factor 5 is defined as the ratio of the moment de- 

veloped at the far edge y= (fig. 9) to the moment at the 
near edge y=~T£ The moment developed at the far edge is 

(A13, 


where w is the deflection of the plate of figure 9, given by 
equation (AlO). 

If the indicated differentiation of equation (AlO) is made 
and the result substituted in equation (A 13), it is found that 



a , 1 a , 0 , 0 a a , 0 . 0 

2 tanh 2 + 2 tan 2“ 2 COth 2 + 2 Cot 2 lf « 

a , a , 0 . 0,a , a 0 0 0 C ° 3 X 

2 ^ 2 + 2 tan 2+2 COth 2“2 Cot 2 

(A 14) 
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The moment at the near edge is, from equation (A8), 

(ill), * =M 0 cos ? (A15) 

'■ 2 A 

* 

By definition, the carry-over factor is 
c _ - (AQ, .a (f tanh f+f tan |)-(f coth f-f cot |) 

(f tanl1 1+ f tan |) +(f coth f-| cot |) 

(A16) 

with the sign of the moment at the far edge changed to 
conform to the sign convention given in the section 
“Definitions." 





Stiffness of a plate with its far edge hinged. — Figure 10 
shows a flat rectangular plate under compression with the 

two edges y= binged to supports. A moment M is 

applied to the edge y =— and the stiffness of the plate is 
defined as 



(A17) 


where (6) b is the rotation of the edge y= — ~ expressed ii 

f- , * 

quarter-radians. The general expression (A5) will again b« 
used to compute 8 and the boundary conditions will be: 


( W ) L =0 
r *2 


(A18 



(A20) 


By use of these boundary conditions, the arbitrary con- 
stants in equation (A5) may be computed, and it is found 
that the deflection surface for the case of figure 10 is 


w= 


Mob 2 


2D{o?+(?) 


sinh 


ay 

T 


cosh nr COS TT sin 


fly 


0y 

T 


sinh ~ cosh ^ 


0 ■ P 

cos ^ sm 2 


rX 

C0S T 


(A21) 


From this deflection surface, the magnitude of the rota- 
tion 8 along the edge y=— ^ is found to be 


(s) r _4 

2 Mb / a s .a B\ 

~W(^T0 r )\ a tanh 2+0 tan 2 +a COth 2~ ? COt V 

(A22) 

where 8 is expressed in quarter-radians. Upon substitution 
of this expression for 8 in equation (Al7), it is found that 

2- p (IKI)' (A23) 

T> a , , a ,0 f 01 a , a 0 0 ' 

2 tani 2+2 tan 2 + 2 COth 2 2 Cot 2 

According to the boundary condition given in equation 

(A20), there is no moment at the edge y—^ Hence, the 

carry-over factor C 11 with the far edge hinged is zero. 

. Stiffness of a plate with far edge free. — Figure 11 shows 
a fiat rectangular plate under compression with one edge 
y=b free and a moment M applied to the parallel edge 
y= 0. The stiffness of the plate is defined as 

<A24) 

where (0)„_ 0 is the rotation of the plate along the edge y — 0 
and is expressed in quarter-radians. 

The general expression (A 5) is used to compute the 
rotation 8. The boundary conditions for a plate with far 
edge free are: 


* 

1 

D 

II 

o 

(A25) 


(A26) 


(A27) 

+ 

1 

V 

II 

o 



=M=Mo cos ^ 
b A 


(A19) 



+ (2-m) 



(A28) 
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Upon determination of the arbitrary constants in equation (A5) by use of these boundary conditions, the deflection 
surface for the case of figure 1 1 is found to be 


where 


>(«*+< s 2 ) 


J^cos cosh sinh — 

-( f )'-- (0 


7 ia (sinh. a— 3> cosh a)— mft sin ft . fty - ] rx 


mi 3 cos ft 


H nx 

j cos T 


(A29) 


/ftY , /x6V 
n_ U/ +/ A2X/ 


m 2 ft cosh a cos ft— n 2 a sinh_« sin p-\-mnp 


m?p smn a cos ft— wa cosh a sin ft 


From the deflection surface, the rotation along the edge y=0 is found to be 

(a\ , /bw\ [~ 4 Adb ~| 2aPmn+aP(m i +n i ) cosh a cos ft-f sinh a sin P .. . 

Wy-o— ■ LS^+jS 2 ) J m 2 ft sinh"a"cos ft— n* ft cosh a sin ft ' ' 

whore 8 is expressed in quarter-radians. Upon substitution of this expression for 8 in equation (A24), the stiffness is found to bo 


^“=T[(l)’+(f 


2m "f f 

cosh 2 H cos 2 ^ 


— m 1 ^1— tan 2 0 tanh j n 2 ^1-ftanh 2 tan ^ 

"(I) (I) ( m2 + n 0 ( ! + tan h 2 f ) (l - tan* |)+4 [m 2 (f )- n 2 (§)* tanh f tan | 


(A3 1 ) 


The trigonometric and hyperbolic functions have been converted to~the half angle in order that the same functions can bo 
used as in the calculation of the other stiffnesses. 


According to the boundaiy condition of equation (A27), 
there is no moment along the edge y—b. The carry-over 
factor G 111 is thus zero for the far edge free. 


Stiffness of a plate with equal and opposite moments 
applied along the unloaded edges. — Figure 12 shows a flat 
rectangular plate under compression, with equal and opposite 



Figure II.— Flats with moment applied at near edge, far edge flee. 


plate is defined as 

s ‘ v= (t),~4 (A32) 

where (6) ^ is the rotation along the edge y —~\ expressed 

*”~2 £ 
in quarter-radians. 


tar 



Figure 12.— Plate with moment applied at near edge, equal and opposite moment 

at far edge. 



PRINCIPLES OF MOMENT DISTRIBUTION APPLIED TO STABILITY OF STRUCTURES COMPOSED OF BARS OR PLATES 69 


The boundary conditions for this case are: 

(w) , 6=0 (A33) 

r=± 2 

V (w +I &)^— W =-‘ V - C ° S T < A34 > 

According to the sign convention of the appendix, the mo- 
ments at the two edges have the same sign although they 
act in opposite directions. By means of these boundary 
conditions, the arbitrary constants in equation (A5) may 
be computed, and the deflection surface for the case of 
figure 12 is found to be 


plate of figure 8 is reversed, the plate will be under tension 
and equation (Al) -will become 


bhc , n dh/? , bho AJc cAw n 
bAby* 1 by* l r ad~ u 


(A38) 


The formal solution of this equation is precisely the same as 
equation (A5), except that the parameters a and 8 are now 
defined by 



(A3 9) 


(A40) 


V}—= 


Mob 2 


D&+0*) 


f" 8y . 

COS -y cosh -y 

cos ^ cosh ^ 


cos Y (A35) 


From this deflection surface, the rotation 8 along the edge 
y——^ is found to be 


Because the stiffnesses S, S 11 , S IU , and5 IV , and the carry- 
over factor C, as calculated for a plate under compression, 
are based directly upon equation (A5), it follows that the 
expression derived for each one of these quantities is still 
correct when the plate is under tension, provided a and 8 
are now given by equations (A39) and (A40). 

The new expressions for a and 8 are complex and may be 
written in the form 


2 Mb 1 

D 

1 

(!) 

1+1 

(I) 

’J 


(A36) 

which is expressed in quarter-radians. Substitution of this 
expression for 8 in equation (A32) gives the stiffness of the 
plate, 

r D (f) + (f) 5“ 




'26tt, , a 8 t fl T+C? 

9 tanh 9+9 tan 9 


CA37) 


Because the moment at is equal and opposite to that 
at y=—^> the carry-over factor C Iv is —1 in accordance with 
the sign convention given in the section “Definitions.” 

PLATE UNDEfi TENSION 

If the direction of the applied longitudinal force on the 


|=^+iB (A41) 

| =B+iA (A42) 

where 

A ~\ yfk V +i+ l (A43) 

b -\ VSVV(D ,+i -r <A44) 

The expressions (A41) and (A42) for a and dare substituted 
into equations (Al 2 ) for S, equation (A23) for S u , equation 
(A31) for S IU , equation (A37) for S 1V , and equation (Al 6 ) 
for C. The results of the substitution show that, for a plate 
in tension, 

0 D AB A sin 45—5 sinh 4 A 
^~b 2 A 1 sin 3 25-5* sinh 5 2+ 

nr ..Pa i p (cosh2A+ cos 25) (sinh A+ sin 2 5) 
b 5 sinh 4 A— A sin 4 B 


(A45) 

CA46) 


am d A(m 2 —l5A 1 B lJ r%TnB T ) sin 45— 5(m 2 — IBAAB 2 — 8m+ J ) sinh 4+ , , 

S == F A£ T&( : m-±A*)*-A*(m+4B i ) i (A4/ ) 

■j — (+ 2 +B 2 ) (m 2 — l&A^BP) (cosh 2 2 A cos 2 25+sinh 2 2 A sin 2 2 B) 

| +[A 2 (m+45 2 ) 2 — B^m— AA^^Ksinh 2 2A cos 2 25+ cosh 2 2 A sin 2 25) 


where 


c „ v D . D cosh 2++cos 25 
5 = 5" jA °5"sinh 2l4 + + sin 25 

„ A cosh 2 A sin 2 5—5 sinh 2 A cos 25 
~B cosh 2A sinh 2 A— A sin 25 cos 25 

m=4(A 2 -5 2 )- M (Yy 


(A48) 

(A49) 

(A50) 


These formulas permit tables of stiffnesses and carry-over 
factors to be prepared for a plate in tension similar to the 
tables of reference 5 for a plate in compression. Such tables, 
however, have not been prepared, and in lieu of them, 
formulas (A45) to (A49) may be used directly if the need 
should arise. 


S4310T — 50- 


■a 
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Members be and it 


(Ibtoin.) 0b/sqln.) 


TABLE I 

RESULTS FOR SOLUTION" OF BAR PROBLEM > 

[For member tf, P=0, P/A -0, E- 28X10» Ib/sq In., 0, and fin,/- 3-397X10* lb-ln.] 


Sic S*ic (Pic p, SP^dc 

(Ib-fn.) (lb-in.) 0b-in.) 


82,700 j tt 10 


1 Table Is from reference 8, table in. 



5. 45X10* 
a 24 

6. 07 
4. 99 
4.94 

4.88 
4.92 
4.99 

aia 

8.34 

6.84 

8.88 
8.33 
a 91 
7.71 
8.89 
0.98 


6.07X10* 

6.01 

4.95 

4.90 

4.86 

4.80 

4.76 

4.72 

4.G8 

4.85 

4.02 

4.68 

4.66 

4.63 

4.60 

4.48 


-2.49X10* 
-3.41 
-4.66 
-a 44 
- 10.01 
-16.5# 

— 51.71 
49.12 
16.88 
10 26 
7.81 

8.73 
4.67 

3.73 
3.098 
2.633 
X 239 



TABLE II 

RESULTS FOR SOLUTION OF PLATE PROBLEM 


JP lr ffl T r /tr\t/h.\ S^r 

(3/6) w\tr J \EwJ (D/b/r 


(kw) 0 _o 

(step #) 


F r ~ 13 ' 


0.8 

3.2 ... 

L 6820 

08410 

0 7555 

1.6965 

.0 

8.6 

L4838 

.7419 

.5994 

1.3413 

L0 

. 4.0 

i.2370 

.6185 

.3729 

.9914 

L2 

4.8 

.4937 

.2468 

—.6729 

-.4261 




0.7 

I 

2.8 

0.5907 

a 2954 

-0.0621 

02333 

.8 

3.2 

.4221 

.2111 

-.3040 

-.0929 

.9 

3.6 

.2275 

.U37 

-.6548 

-.5411 

1.0 

4.0 

0 

0 

-L2117 

-1.2117 


FT*' 0 ' £- 1.8 


0.7 

2.3 

04373 

0.21S6 

-a 1640 

0.0516 

.8 

3.2 

.3568 

.1784 

-.2868 

-.1084 

.9 

3.6 

.2706 

.1353 

-.4351 

-.2998 

L0 

4.0 

.1781 

.0890 

—.6179 

-.5289 


&-**• FT™ 


-01266 01058 

-,1799 . 0205 9 

-.2501 -.0726 S - Z * B 

-.3285 -.1747 






























